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Abstract. Specifying boundary conditions continues to be a challenge in nu- 
merical relativity in order to obtain a long time convergent numerical simu- 
lation of Einstein's equations in domains with artificial boundaries. In this 
paper, we address this problem for the Einstein-Christoffel (EC) symmetric 
hyperbolic formulation of Einstein's equations linearized around flat space- 
time. First, we prescribe simple boundary conditions that make the problem 
well posed and preserve the constraints. Next, we indicate boundary conditions 
for a system that extends the linearized EC system by including the momen- 
tum constraints and whose solution solves Einstein's equations in a bounded 
domain. 



1. Introduction 

In the Arnowitt-Deser-Misner or ADM decomposition, Einstein's equations split 
into a set of evolution equations and a set of constraint equations (see Section [5]) , 
and what one does to construct a solution consists of first specifying the initial data 
that satisfies the constraints and then applying the evolution equations to compute 
the solution for later times. The problem of well-posedness in the analytic sense 
has been intensely studied, with the result that there is a great deal of choice of 
formulations available for analytic studies (see [15], [l^, [29], [6], [32], [1], [I], [8], 
[I], [n], [13], [H], [20], [21], [23], [27], [40], [H], among others). However, in numer- 
ical relativity, one usually solves the Einstein equations in a bounded domain (cubic 
boxes are commonly used) and the question that arises is what boundary conditions 
to provide at the artificial boundary. In general, most numerical approaches have 
been made using carefully chosen initial data that satisfies the constraints. On the 
other hand, finding appropriate boundary conditions that lead to well-posedness 
and consistent with constraints is a difficult problem and subject to intense in- 
vestigations in the recent years. In 1998, Stewart |49] has addressed this subject 
within Frittelli-Reula formulation [29) linearized around flat space with unit lapse 
and zero shift in the quarter plane. Both main system and constraints propagate 
as first order strongly hyperbolic systems. This implies that vanishing values of 
the constraints at t = will propagate along characteristics. One wants the values 
of the incoming constraints at the boundary to vanish. However, one can not just 
impose them to vanish along the boundaries since the constraints involve deriva- 
tives of the fields across the boundary, not just the values of the fields themselves. 
If the Laplace-Fourier transforms are used, the linearity of the differential equa- 
tions gives algebraic equations for the transforms of the fields. Stewart deduces 
boundary conditions for the main system in terms of Laplace-Fourier transforms 
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that preserve the constraints by imposing the incoming modes for the system of 
constraints to vanish and translating these conditions in terms of Laplace-Fourier 
transforms of the main system variables. In 1999, a well posed initial-boundary 
value formulation was given by Friedrich and Nagy [22j in terms of a tetrad-based 
Einstein-Bianchi formulation. In view of our work which is to be presented here, 
of particular interest are the more recent investigations regarding special boundary 
conditions that prevent the influx of constraint violating modes into the computa- 
tional domain for various hyperbolic formulations of Einstein's equations (see [2] , 
[31, 0, [H], US], [17], [SI, [S], [ig, [H], [ig, [SH, [Sg, among others). A differ- 
ent approach can be found in [28], where the authors stray away from the general 
trend of seeking to impose the constraints along the boundary. They argue that 
the projection of the Einstein equations along the normal to the boundary yields 
necessary and appropriate boundary conditions for a wide class of equivalent for- 
mulations. The ideas and techniques introduced in |28j are further developed and 
proven to be effective by the same authors in [25] . In principal, they show that the 
projection of the Einstein tensor along the normal to the boundary relates to the 
propagation of the constraints for two representations of Einstein's equations with 
vanishing shift vector, namely, the Arnowitt-Deser-Misner (ADM) formulation [TO] 
and the classical Einstein-Christoffel (EC) formulation |5]. In particular, for the 
EC formulation they retrieve a subclass of the boundary conditions presented in 
[17] and [9|. 

Of course, specifying constraint-preserving boundary conditions for a certain 
formulation of Einstein's equations does not solve entirely the complicated problem 
of numerical relativity. There are other aspects that have to be addressed in order 
to obtain good numerical simulations; for example, the existence of bulk constraint 
violations, in which existing violations are amplified by the evolution equations (see 
[Is] . [19], [36], [45], and references therein). A review of some work done in this 
direction can be found in the introductory section of |34j . Before we end this very 
brief review, it should also be mentioned the work done on boundary conditions for 
Einstein's equations in harmonic coordinates, when Einstein's equations become 
a system of second order hyperbolic equations for the metric components. The 
question of the constraints preservation does not appear here, as it is hidden in the 
gauge choice, i.e., the constraints have to be satisfied only at the initial surface, 
the harmonic gauge guarantees their preservation in time (see [38], [50], [51], and 
references therein). 

In this paper we address the boundary conditions problem for the classical EC 
equations derived in [6], linearized with respect to the flat Minkowski spacetime, 
and with arbitrary lapse density and shift perturbations. This problem has been 
addressed before in [16] in the case of spherically symmetric black-hole spacetimes 
in vacuum or with a minimally coupled scalar field, within the EC formulation of 
Einstein's equations. Here Stewart's idea of imposing the vanishing of the ingo- 
ing constraint modes as boundary conditions is employed once again. Then, the 
radial derivative is eliminated in favour of the time derivative in the expression of 
the ingoing constraints by using the main evolution system. The emerging set of 
boundary conditions depends only on the main variables and their time derivative 
and preserves the constraints. In [17j this technique is refined and employed for 
the generalized EC formulation [32] when linearized around Minkowski spacetime 
with vanishing lapse and shift perturbations on a cubic box. Again, the procedure 
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consists in choosing well-posed boundary conditions to the evolution system for 
the constraint variables and translating them into well-posed boundary conditions 
for the variables of the main evolution system. The scheme proposed in [17 ends 
up giving two sets, called "Dirichlet and Neumann-like," of constraint preserving 
boundary conditions. However, the energy method used in [17 works only for 
symmetric hyperbolic constraint propagation, which forces the parameter rj of the 
generalized EC system to satisfy the condition < 77 < 2. Therefore the analysis 
in [T7] does not cover the case = 4 required for the standard EC formulation 
introduced in [IT. In [5] we announced and presented our results on the boundary 
conditions problem for the standard EC formulation (77 = 4) linearized around the 
Minkowski spacetime with arbitrary lapse density and shift perturbations in the 
Penn State Numerical Relativity Seminar. In essence, we introduced the very same 
sets of boundary conditions that are under scrutiny in this material, i.e., ()18p and 
(P7)) (see [5]). Much of this material appeared also in the thesis of the second author 

m- 

The organization of this paper is as follows: in Section [2] we introduce Ein- 
stein's equations and their ADM equations for vacuum spacetime. In Section [3l 
by densitizing the lapse, linearizing, and defining a set of new variables, we derive 
the linearized EC first order symmetric hyperbolic formulation around flat space- 
time. The equivalence of this formulation with the linearized ADM is proven in 
the Cauchy problem case. In Section [4] we indicate two distinct sets of well-posed 
constraint-preserving boundary conditions for the linearized EC. We prove that the 
linearized EC together with these boundary conditions is equivalent with linearized 
ADM on polyhedral domains. In Section [5] we indicate boundary conditions for an 
extended unconstrained system equivalent to the linearized ADM decomposition. 
In Section [6] we discuss the case of inhomogeneous boundary conditions. We end 
this work with a summary and a discussion of our results in Section [T] For reader's 
convenience, in the appendix we review a classical result on the well-posedness 
of maximal nonnegative boundary conditions for symmetric hyperbolic systems. 

2. Einstein's Equations and the ADM Decomposition 

In general relativity, spacetime is a 4-dimensional manifold M of events endowed 
with a pseudo-Riemannian metric gajs that determines the length of the line ele- 
ment ds"^ = ga(3dx"dx^ . This metric determines curvature on the manifold, and 
Einstein's equations relate the curvature at a point of spacetime to the mass-energy 
there: Ga[3 — SttTq^, where Gap is the Einstein tensor, i.e., the trace-reversed Kicci 
tensor Gap ■— Rap — ■^Rgap, and Tap is the energy-momentum tensor. In what 
follows we will restrict ourselfs to the case of vacuum spacetime, that is Tap — 0. 
Einstein's equations can be viewed as equations for geometries, that is, their so- 
lutions are equivalent classes under spacetime diffeomorphisms of metric tensors. 
To break this diffeomorphisms invariance, Einstein's equations must be first trans- 
formed into a system having a well-posed Cauchy problem. In other words, the 
spacetime is foliated and each slice St is characterized by its intrinsic geometry 7^ 
and extrinsic curvature Kij , which is essentially the "velocity" of jij in the unit 
normal direction to the slice. Subsequent slices are connected via the lapse function 
N and shift vector /3* corresponding to the ADM decomposition [10] (also [54]) of 
the line element 



ds^ = -N^dt^ + 7^J(dx'' + (3'dt){dx^ + P'dt). 



(1) 
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This decomposition allows one to express six of the ten components of Einstein's 
equations in vacuum as a constrained system of evolution equations for the metric 
7y and the extrinsic curvature Kij: 

k,, - N[R,, + {K\)K,j - 2KaK'-j\ + V^Xy + KuV.fi' + KijV,l3' - V,V,7V, 

V'K,, - V,K] = 0. 

(2) 

where we use a dot to denote time differentiation and Vj for the covariant derivative 
associated to ^ij . The spatial Ricci tensor i?y has components given by second order 
spatial differential operators applied to the spatial metric components 7ij. Indices 
are raised and traces taken with respect to the spatial metric 7^ , and paranthesized 
indices are used to denote the symmetric part of a tensor. 

3. Linearized Einstein-Christoffel 

The Einstein-Christoffel or EC formulation |6j is derived from the ADM system 
with a densitized lapse. That is, we replace the lapse N in ^ with a^/g where a 
denotes the lapse density. A trivial solution to this system is Minkowski spacetime 
in Cartesian coordinates, given by jij = Sij, Kij = 0, = 0, a = I. In the 
remainder of the paper we will consider the problem linearized about this solution. 
To derive the linearization, we write 7^ = Sij + cjij, Kij — Kij, (3^ — (3'' , a — 1 + a, 
where the bars indicate perturbations, assumed to be small. If we substitute these 
expressions into ([2]) (with N = a^J^), and ignore terms which are at least quadratic 
in the perturbations and their derivatives, then we obtain a linear system for the 
perturbations. Dropping the bars, the system is 

= -2if,, +2%/3j), (3) 
= d^d{j9i)i - \d^dig,j - d,djg\ - d.djU, (4) 
C:=9^(5'.gy-9,.g|) = 0, (5) 
Cj := d^Kij - djK\ = 0, (6) 
where we use a dot to denote time differentiation. 

Remark. For the linear system the effect of densitizing the lapse is to change the 
coefficient of the term didjg\ in Had we not densitized, the coefficient would 
have been —1/2 instead of —1, and the derivation of the linearized EC formulation 
below would not be possible. 

The usual approach to solving the system ([3|)-([6]) is to begin with initial data 
gij{Q) and i^ij(O) defined on and satisfying the constraint equations ([5|), ([6]), 
and to define gij and Kij for i > via the Cauchy problem for the evolution 
equations ([2]), (H). It can be easily shown that the constraints are then satisfied for 
all times. Indeed, if we apply the Hamiltonian constraint operator defined in ([5]) to 
the evolution equation ([3]) and apply the momentum constraint operator defined in 
^ to the evolution equation ([5]), we obtain the first order symmetric hyperbolic 
system 

C = -2dW„ C,^~\dfi. 
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Thus if C and Cj vanish at < = 0, they vanish for all time. 

The linearized EC formulation provides an alternate approach to obtaining a 
solution of ([2])^® with the given initial data, based on solving a system with 
better hyperbohcity properties. If gij, Kij solve ©-(ini), define 

fki] = ]^[dkgij - [d'-gu - d^g'i)Sjk - {d^gij - djg\)5,k]- (7) 
Then —d^ fkij coincides with the first three terms of the right-hand side of (|3]), so 

= -^''^k^J - d,d,a. (8) 

Differentiating ([7]) in time, substituting ([3]), and using the constraint equation ([5]), 
we obtain 

fkij = —dkKij + Lkij, (9) 

where 

Lktj = dkd^ifSjf ~ d^d[i(3i]Sjk - d^d[iPj]6,k (10) 

The evolution equations ([8]) and ([9]) for Kij and fkij, together with the evolution 
equation ([3]) for gij, form the linearized EC system. As initial data for this system 
we use the given initial values of gij and Kij and derive the initial values for fkij 
from those of gij based on 

fkrjiO) - l{^kg^m - [d'gim - d,g\mS,k ~ [d'gijiO) - d^gHoMk}. (11) 

In this paper we study the preservation of constraints by the linearized EC system 
and the closely related question of the equivalence of that system and the linearized 
ADM system. Our main interest is in the case when the spatial domain is bounded 
and appropriate boundary conditions are imposed, but first we consider the result 
for the pure Cauchy problem in the remainder of this section. 

Suppose that Kij and fkij satisfy the evolution equations ([8]) and ^ (which 
decouple from ([3])). If Kij satisfies the momentum constraint ([6]) for all time, then 
from ([5]) we obtain a constraint which must be satisfied by fkij '■ 

d\d'fki,-d,fki') = o. (12) 

Note that ([7]) is another constraint that must be satisfied for all time. The following 
theorem shows that the constraints (H]) , ([7]) , and ([12]) are preserved by the linearized 
EC evolution. 

Theorem 1. Let initial data gij{0) and Kij{0) be given satisfying the constraints 
^ and respectively, and fkij{0) be defined by (llip . Then the unique solution 
of the evolution equations and ^ satisfies ([S]), ([7]), and (jl2p for all time. 

Proof. First we show that the initial data fkij{0) defined in pTjl satisfies the con- 
straint Applying the constraint operator in to (fTT|) we find 

d'id'fkij - d,fk/W) = \d,{d'd'gki - d'^dkglm - ld,C{Q) = 0. (by ©) 

It is immediate from the evolution equations that each component Kij satisfies the 
inhomogeneous wave equation 

K,, = d'^dkK,, - d^'Luj - d^dja. 
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Applying the momentum constraint operator defined in ([6|) , we see that each com- 
ponent Cj satisfies the homogeneous wave equation 

Cj = d'^dkCj. (13) 

Now Cj = at the initial time by assumption, so if we can show that Cj = at 
the initial time, we can conclude that Cj vanishes for all time. But, from ([8]) and 
the definition of Cj , 

C,^~dHd'hi,~d,hi'l (14) 

which we just proved that vanishes at the initial time. Thus we have shown Cj 
vanishes for all time, i.e., fH]) holds. In view of (1141) . (fT2|) holds as well. From ^ 
and ([U we have 

fkij = ^dkgij - d'd[if]i]6jk - d'd[iP^]6ik. 
Applying the momentum constraint operator to ([3]) and using ([6|), it follows that 

so fkij - [dkgij - [d^gii - dig\)5kj - (d'gij - djgl)Ski]/2 does not depend on time. 
From ((ni), we have □ 

In view of this theorem it is straightforward to establish the key result that for 
given initial data satisfying the constraints, the unique solution of the linearized 
EC evolution equations satisfies the linearized ADM system, and so the linearized 
ADM system and the linearized EC system are equivalent. 

Theorem 2. Suppose that initial data gij(0) and Kij(0) are given satisfying the 
Hamiltonian constraint ^ and momentum constraint respectively, and that 
initial data fkij (0) is defined by (jlll) . Then the unique solution of the linearized 
EC evolution equations ^ satisfies the linearized ADM system 

Proof. From Theorem [1] we know that Cj = for all time, i.e., (O holds. Then 
from ^ and ^ we see that C — ~2d^Cj — 0, and, since C vanishes at initial time 
by assumption, C vanishes for all time, i.e., (O holds as well. 

It remains to verify Q. From Theorem [H we also have (O. Substituting ([7]) in 
dl]) gives (HJ, as desired. □ 

4. Maximal Nonnegative Constraint Preserving Boundary Conditions 

In this section of the paper, we provide maximal nonnegative boundary condi- 
tions for the linearized EC system which are constraint-preserving in the sense that 
the analogue of Theorem[l]is true for the initial-boundary value problem. This will 
then imply the analogue of Theorem [2l We assume that f2 is a polyhedral domain. 

Consider an arbitrary face of dfl and let denote its exterior unit normal. 
Denote by m* and two additional vectors which together form an orthonormal 
basis. The projection operator orthogonal to is then given by := mim^ + liP 
(and does not depend on the particular choice of these tangential vectors). Note 
that 

Si=n,n^+T^, T^r^=r^. (15) 

Consequently, 

viw'' — n-'vjUiW^ + rlvjrlw^ for all u/, vJ" . (16) 
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First we consider the following boundary conditions on the face: 

n'm'K^j = n'-PK^j = n'^n'n^fMj = n'^rn'm^ f^j = nf^rP fkio = n'^m^^ fkzj = 0. 

(17) 

These can be written as well: 

nW^''K,j = 0, n''nVfk^J = 0, tiV^'t^"^,, = 0, (18) 

and so do not depend on the choice of basis for the tangent space. We begin by 
showing that these boundary conditions are maximal nonnegative for the hyper- 
bolic system ([3]), ([8]), and ([9]), and so, according to the classical theory of [ST and 
[35j (also [30], [33], [37], [39], [46], [47], among others), the initial-boundary value 
problem is well-posed. For convenience, in Appendix A we recall the definition and 
a classical result due to Ranch [39] on well-posedness of maximal nonnegative 
boundary conditions. 

Let V denote the vector space of triplets of constant tensors {gij,Kij, fkij) all 
three symmetric with respect to the indices i and j. Thus dimV = 30. The 
boundary operator An associated to the evolution equations ([3]) , ([5]) , and ^ is the 
symmetric linear operator V V given by 

fjij = 0, K^j = n'^fkij, fki] = nkK.j. (19) 

A subspace iV of is called nonnegative for An if 

g^Jg'' + K,, + h^, f"'' > (20) 

whenever {gij,Kij, fkij) G N and (g^, Kij, fkij) is defined by (|19p . The subspace is 
maximal nonnegative if also no larger subspace has this property. Since An has six 
positive, 18 zero, and six negative eigenvalues, a nonnegative subspace is maximal 
nonnegative if and only if it has dimension 24. Our claim is that the subspace N 
defined by (|17p is maximal nonnegative. The dimension is clearly 24. In view of 
p^ . the verification of (|20p reduces to showing that fkijK'^-' > whenever pTj) 
holds. In fact, fkijK^-' — 0, that is, n^fkij and Kij are orthogonal (when (fT7|) 
holds). To see this, we use orthogonal expansions of each based on the normal and 
tangential components: 

Kij = n'-mrrrijKin^ + n^n.T^Kim -t- rliiJ^njKhn + r^r^if/™, (21) 

n'^fktj = n^rnn'^njn^fkim ^ n^riiT^'n'' fkim + rlrf^n^n'' fki„, + rlrj'n^fkim. (22) 

In view of the boundary conditions (in the form (118^ ). the two inner terms on the 
right-hand side of (|2T]) and the two outer terms on the right-hand side of ([22]) 
vanish, and so the orthogonality is evident. 

Next we show that the boundary conditions are constraint-preserving. This is 
based on the following lemma. 

Lemma 3. Suppose that a and (3^ vanish. Let gij, Kij, and fkij be a solution to 
the homogeneous hyperbolic system ([3]), ([8]), and ^ and suppose that the boundary 
conditions (I17p are satisfied on some face of dQ. Let Cj be defined by ([6]). Then 

CjU^diC^ = (23) 

on the face. 

Proof. In fact we shall show that n^Cj = (so also n^Cj — 0) and r^n^diC^ = 0, 
which, by (fTB]) implies (p3| . First note that 

C, = (<5f <5^'= - <5|,5"")9feX™ = (5f nV'' + S^'t'" - Sp'"')dkK,m, 
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where we have used the first identity in (|15|) . Contracting with gives 

= -n n n jkim + t Ti n duKim + n d fkim 

where now we have used the equation ^ (with j3i — 0) for the first and last 
term and the second identity in for the middle term. From the boundary 
conditions we know that n^n^n^ fkim = 0, and so the first term on the right-hand 
side vanishes. Similarly, we know that r'^^n'^Kim = on the boundary face, and 
so the second term vanishes as well (since the differential operator rfdk is purely 
tangential). Finally, n''S^™'fkim = n^in^n^ + ri™ + m^rn^)fkim = 0, and so the 
third term vanishes. We have established that n^Cj = holds on the face. 
To show that r^nJ'diC^ = on the face, we start with the identity 

Similarly 
Therefore, 

= {rP'^n'S''^ - tP'^u't^^ + r^'^T'^n^ - rP^n^n'ri"' - TP^T'"'n^)dkdiK^„,. 

For the last three terms, we again use ([9]) to replace diKim with —fum and argue as 
before to see that these terms vanish. For the first term we notice that S^^dkdiKim = 
d^dkKim = Kim (from ^ and (O with vanishing a and Since r'P'^n'^Kim 
vanishes on the boundary, this term vanishes. Finally we recognize that the second 
term is the tangential Laplacian, r^^dkdi applied to the quantity li^TP^Kim, which 
vanishes. This concludes the proof of (|23|) . □ 

The next theorem asserts that the boundary conditions are constraint-preserving. 

Theorem 4. LetQ be a polyhedral domain. Given gijiO) and Kijifi) onD, satisfying 
the constraints ([5]) and respectively, and fkiji^) defined by (jlip . define gij, 
Kij, and fkij for positive time by the evolution equations ([5]), and ^ and the 
boundary conditions (jl7p . Then the constraints ([5]), ([7]), and (|12p are satisfied for 
all time. 

Proof. Exactly as for Theorem [1] we find that Cj satisfies the wave equation 
and both Cj and Cj vanish at the initial time; these facts are unrelated to the 
boundary conditions. Define the usual energy 

E{t) ^ \ I {CjC^ + d'CjdiC^)dx. 

Clearly E{0) ~ 0. From (fT^)) and integration by parts 

[ CjU^diC'da. (24) 

Therefore, if a = and /?' = 0, we can invoke Lemma [31 and conclude that E is 
constant in time. Hence E vanishes identically. Thus Cj is constant, and, since it 
vanishes at time 0, it vanishes for all time. By (fT4|) . the constraints (fT2|) are also 
satisfied for all time. This establishes that the constraints ([6|) and (fT2| hold under 
the additional assumption that a and /3' vanish. 
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To extend to the case of general a and /3' we use Duliamel's principle. Let 
S{t) denote the solution operator associated to the homogeneous boundary value 
problem. That is, given functions /lij(O), Ky(0), (j)kij{^) on fi, define 

5(t)(/ly(0),K,,(0),<^fe„-(0)) = {h,j{t),K^j{t),d^kr3{t)), 

where hij , Hij , (j)kij is the solution to the homogeneous evolution equations 

hij — 2k^j, Kij — d ' (pkij : 4*kij — ^k^ij : 

satisfying the boundary conditions and assuming the given initial values. Then 
Duhamel's principle represents the solution Qij, Kij, fkij of the inhomogeneous 
initial-boundary value problem ([5]), (O, PT|) as 

+ f S{t~s){2d^,l3^),-d,d,a{s),Lk^j{s))ds. (25) 

JQ 

Now it is easy to check that the Hamiltonian constraint (O is satisfied when gij is 
replaced by 29(^/3^) (for any smooth vector function /3'), the momentum constraint 
([6]) is satisfied when Kij is replaced by —didja{s) (for any smooth function a), and 
the constraint (fT2|) is satisfied when fkij is replaced by Lkij (s) defined by (fTO|) (for 
any smooth vector function /3*). Hence the integrand in (|25p satisfies the constraints 
by the result for the homogeneous case, as does the first term on the right-hand 
side, and thus the constraints ^ and (fT2|) are indeed satisfied by Kij and fkij, 
respectively. 

The proof of the fact that the constraints ([7]) are satisfied for all time follows 
exactly as in Theorem [T] 

Note that the boundary conditions (|f 7p play a crucial role in proving that the 
momentum constraints ^ are preserved for all time; the preservation of the con- 
straints ([7]) and (|f 2p being a consequence of this fact. □ 

The analogue of Theorem [2] for the initial-boundary value problem follows from 
the preceding theorem exactly as before. 

Theorem 5. Let Q be a polyhedral domain. Suppose that initial data gtjiO) and 
Kij{0) are given satisfying the Hamiltonian constraint ([5]) and momentum con- 
straint Q , respectively, and that initial data fkij (0) is defined by (|f f |) . Then the 
unique solution of the linearized EC initial-boundary value problem ([8]), 
together with the boundary conditions (|f 7p satisfies the linearized ADM system ([5])- 
^ in n. 

We close this section by noting a second set of boundary conditions which are 
maximal nonnegative and constraint-preserving. These are 

n'n^Kij = m'm^Kij = PPK.j = m'PKij = n'^n'm^ fk,j = n^n'P fk^J = 0, (26) 

or, equivalently, 

nVif,j=0, r^'r^^if,, = 0, nVr^Vfe..=0. (27) 

Now when we make an orthogonal expansion as in (pij) . (P^ . the outer terms on 
the right-hand side of the first equation and the inner terms on the right-hand side 
of the second equation vanish (it was the reverse before), so we again have the 
necessary orthogonality to demonstrate that the boundary conditions are maximal 
nonnegative. Similarly, to prove the analogue of Lemma [3l for these boundary 
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conditions we show that the tangential component of Cj vanishes and the normal 
component of n^diC^ vanishes (it was the reverse before). Otherwise the analysis 
is essentially the same as for the boundary conditions pT|) . 

5. Extended EC System 

In this section we indicate an extended initial boundary value problem whose 
solution solves the linearized ADM system Q-® in This approach could 
present advantages from the numerical point of view since the momentum constraint 
is "built-in," and so controlled for all time. The new system consists of (jS)), ([9]), 
and two new sets of equations corresponding to ([8|) 

= -d'^fk^j + ^{^^PJ + djPi) - d'^pkS^j - didja, (28) 

and to a new three dimensional vector field pi defined by 

p,=d'Ku-d,Kl. (29) 

Observe that the additional terms that appear on the right-hand side of (|28|) com- 
pared with ([8]) are nothing but the negative components of the formal adjoint of 
the momentum constraint operator applied to pi. 

Let V be the vector space of quadruples of constant tensors {gijTKij, fkij^Pk) 
symmetric with respect to the indices i and j. Thus diml^ = 33. The boundary 
operator An ■ V ^ V in this case is given by 

gtj = 0, k^j = n''fkij-^{niPj+njPi)+n''pkSij, fkij = UkKij, pi = -n^ Kii+n,,K\. 

(30) 

The boundary operator A„ associated to the evolution equations ([3]), (I28p . ([H]), and 
(|29p has six positive, 21 zero, and six negative eigenvalues. Therefore, a nonnegative 
subspace is maximal nonnegative if and only if it has dimension 27. We claim that 
the following boundary conditions are maximal nonnegative for ([3]), (|28|) . (I9|), and 

n'^m'K^ = nH^Kij =n^n'^n^ fkij ^ n^irrem^ fktj +Pk) = 
n\VVfk^J + Pk) = n'^m'F/fey = 0. 
These can be written as well: 

nV^'^X,, =0, n'^nV/fe,, =0, n'=(T'V^™/fc,j + r^^) = 0, (32) 

and so do not depend on the choice of basis for the tangent space. 

Let us prove the claim that the subspace N defined by (|3ip is maximal non- 
negative. Obviously, dimA'^ = 27. Hence, it remains to be proven that N is also 
nonnegative. In view of ([30]), the verification of non- negativity of N reduces to 
showing that 

n^ fk^JK'^ - n'p>K^j + n^pkK\ > (33) 
whenever ([3T|) holds. In fact, we can prove that the left-hand side of (|33p vanishes 
pending PT|) holds. From the boundary conditions (in the form ([21])) and the 
orthogonal expansions (|2ip and ([^^ of Kij and fkij , respectively, the first term on 
the right-hand side of ([55]) reduces to n'^r^'-T^^"- fkijKi^ ~ —n'^PkT^"^Kijn- Then, 
combining the first and third terms of the left-hand side of ([33]) gives —n^pkr'^-' Kij + 
n^PkS^'' Kij = n^pkri^n^ Kij. Finally, by using the orthogonal decomposition p' = 
n^pkn^ + T^-'pk and the first part of the boundary conditions ([32]) the second term 
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of the left-hand side of ([33|) is —n'^pkn^n^ Kij —pkn^r^^ Kij ~ ~n''pkn^n^ Kij , which 
is precisely the negative sum of the first and third terms of the left-hand side of 
([55)1 . This concludes the proof of ([55)1 . 

Theorem 6. Let D, be a polyhedral domain. Suppose that the initial data gij(0) 
and Kij (0) are given satisfying the Hamiltonian ([5]) and momentum constraints 
([6]), respectively, fkijiO) is defined by ([7]), and pi{0) = 0. Then the unique solution 
{gij , Kij , fkij ,Pi) of the initial boundary value problem ([3]), (j28p . ([9]), and (j29p . 
together with the boundary conditions (I31|) . satisfies the properties pi — for all 
time, and {gij,Kij) solves the linearized ADM system in f2. 

Proof. Observe that the solution of the initial boundary value problem (|3]), ([8|), ([9]), 
and (|17p (boundary conditions), together with pi — for all time, is the unique 
solution of the initial boundary value problem (pS)) . and together 
with the boundary conditions ([3T|) . The conclusion follows from Theorem [5j □ 

We close by indicating a second set of maximal nonnegative boundary conditions 
(corresponding to ^) for and ^ for which Theorem [6] holds as 

well. These are 

k ^ , k k .„■/ ,fe n (34) 

n'^n'^mP fk.,j - m^pk = n'^n'-P fk^J - tpk = 0, 

or, equivalently, 

nVif„- = 0, r^V^"if„=0, nW/fc^j-rV^O. (35) 

6. Inhomogeneous Boundary Conditions 

In this section we provide a formal method of constructing well-posed constraint- 
preserving inhomogeneous boundary conditions for ([3]), ([5]), and ^ corresponding 
to the two sets of boundary conditions p?]) and (pS)) , respectively. The first set of 
inhomogeneous boundary conditions corresponds to pT)) and can be written in the 
following form 

n'm^k.j = n'PK.j = n'^n'n^fktj = n'^m'm^ f^j = n^PP fkij = n'^m'P fkij = 0, 

(36) 

where Kij = Kij - Kij, fkij = fkij - Fkij, with Kij and given in Vl for aU time 
and satisfying the constraints ([6]) and ^2]) . respectively. 

The analogue of Theorem H] for the inhomogeneous boundary conditions ([36]) is 
true. 

Theorem 7. LetQ be a polyhedral domain. Given gij{Q) and Kij (0) onQ satisfying 
the constraints ([5]) and ([6|), respectively, and fkij{0) defined by (jlip . define gij , 
Kij, and fkij for positive time by the evolution equations ([5]), and ^ and the 
boundary conditions (j36p . Then the constraints ([7]), and (|12p are satisfied for 
all time. 

Proof. Observe that gij, Kij, and /feij satisfy and ([51) with the forcing 

terms replaced by 29(j/3j), —didja — d^Fkij — kij, and L^j^j — dkKij — Fkij, respec- 
tively. Exactly as in Theorem IH it follows that Kij and fkij satisfy ([6]) and p2| . 
respectively, for all time. Thus, Kij and fkij satisfy ([6]) and (fT2|) . respectively, for 
all time. Finally, same arguments as in Theorem [1] show that the constraints ^ 
are also preserved through evolution for all time. □ 



12 



DOUGLAS N. ARNOLD AND NICOLAE TARFULEA 



The analogue of Theorem [5] for the case of the inhomogeneous boundary condi- 
tions ([36|) follows from the preceding theorem by using the same arguments as in 
the proof of Theorem [21 

Theorem 8. Let Q be a polyhedral domain. Suppose that initial data gij (0) and 
Kij (0) are given satisfying the Hamiltonian constraint ^ and momentum con- 
straint ([6|) , respectively, and that initial data fkij (0) is defined by pip . Then the 
unique solution of the linearized EC initial-boundary value problem 
together with the inhomogeneous boundary conditions (p6)) satisfies the linearized 
ADM system ©-([e]) in 

Note that there is a second set of inhomogeneous boundary conditions corre- 
sponding to (j26p for which Theorem [7] and Theorem [5] remain valid. These are 

n'n^Kij = m'm^kij = PPK.j = m'PKij = n'^n'm^ fk,j = n'^n'Pfk.j = 0, (37) 

where again K^j = K^j - k.^, fkij = fkij - Fkij, with and Fkij given and 
satisfying the constraints ([6]) and p2)) . respectively. 

Similar considerations can be made for the extended system introduced in the 
previous section. There are two sets of inhomogeneous boundary conditions for 
which the extended system produces solutions of the linearized ADM system ([3])- 
([S]) on a polyhedral domain fi. These are 

n'm^k.j = n'PK.j = n'^n'n^ fktj = n'' {m'm^ fktj + Pk) = n^{PPfkij + Pk) = 

n^m'PfMj = 

(38) 

and 

n^n^ Kij — m^m^ Kij = m^P Kij — PP Kij = n'^n^m^ fkij — m'^pk — 

n'^nH^fuj - I'^Pk = 0, 

where Kij and fkij are defined as before. 

The next theorem is an extension of Theorem [S] to the case of inhomogeneous 
boundary conditions. 

Theorem 9. Let Q be a polyhedral domain. Suppose that the initial data gij{0) 
and Kij{0) are given satisfying the Hamiltonian ^ and momentum constraints 
respectively, fkijiO) is defined by (O, and pi{0) = 0. Then the unique solution 
{gij , Kij , fkij , Pi) of the initial boundary value problem ([3]), (|28p . ([9|), and (|29p. 
together with the inhomogeneous boundary conditions (|38p (or (j39p ). satisfies the 
properties Pi = for all time, and (gij, Kij) solves the linearized ADM system 

©-® m n. 

Proof. Note that the solution of the initial boundary value problem Q, ([9]), 
and ([36]) (or ([37| . respectively), together with pi = for all time, is the unique 
solution of the initial boundary value problem ([3]), (|28|) . ([9]), and ([29|l . together 
with the boundary conditions ((SS]) (or (j39p . respectively). The conclusion follows 
from Theorem [51 □ 

7. Concluding Remarks 

We have studied the boundary conditions problem for the standard EC formula- 
tion of Einstein's equations linearized about the Minkowski spacetime. In SectionHl 
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we indicate two sets of maximal nomiegative boundary conditions ()17|) and ()26|1 . 
respectively, which are consistent with the constraints. These boundary conditions 
were announced in 9 and overlap with the boundary conditions found in ITJ for 
the generalized EC formulation for < i] < 2 with vanishing shift and lapse density 
perturbations. However, the energy method of [17] works only for the general- 
ized EC formulation with < ry < 2; the standard EC formulation corresponds 
to = 4. Moreover, we prove that our boundary conditions are well-posed and 
consistent with the constraints in the more general case of arbitrary shift and lapse 
density perturbations by using a new argument involving the Duhamel's principle. 
Also, our approach emphasizes the relation between the ADM formulation and the 
constrained evolution of the EC system in the linearized context. Besides, our 
method is simpler, yet effective, and seems to be easily transferable to other formu- 
lations and/or other background spacetimes. In fact, other Einstein's hyperbolic 
formulations, e.g., Alekseenko- Arnold [1], are analyzed in 52J by using the same 
method. A subclass of the boundary conditions presented in this paper and intro- 
duced previously in |9j has been pointed out by Frittelli and Gomez in |25j (in the 
case of vanishing shift vector) as an example to their Einstein boundary conditions, 
that is, the vanishing of the projection of Einstein's tensor along the normal to the 
boundary. 

One of the main results in this paper is the construction of an extended symmetric 
hyperbolic system which incorporates the momentum constrains as main variables. 
For this extended system, we construct two sets of maximal nonnegative boundary 
conditions and establish its relationship with the linearized ADM formulation. Such 
a construction could serve as a model of how to control the bulk constraint violations 
by making the constraints part of the main evolution system, and so keeping them 
under control for all time. To the best of our knowledge, this is a new approach 
regarding the bulk constraint violations control. 

We also make some considerations about how inhomogeneous boundary condi- 
tions consistent with the constraints could be constructed. 

In some places, our methods of proof interfere with the techniques used in |16| 
and |17| . e.g., using the trading of normal derivatives for tangential and temporal 
ones and the use of the energy method to prove that the constraints are preserved. 
We apply these techniques to the slightly more general case of polyhedral domains 
(as opposed to cubic boxes) and in a more systematic way. This could be of po- 
tential interest to the case of curved boundary domains, for which the derivative 
components trading techniques introduce new terms related to the geometry of the 
boundaries (see '52], Section 4.2, for the analysis of a model problem similar to the 
linearized EC formulation on curved domains). It is also expected that these or 
similar techniques will be useful in the nonlinear case. For the interested reader, 
we point out the work done in [34], where the authors construct new boundary 
conditions for the nonlinear KST form [32] of the Einstein equations (which in- 
cludes the EC formulation). Their boundary conditions are designed to prevent 
the influx of constraint violations and physical gravitational waves into the compu- 
tational domain. However, as specified in [3l|, there is no rigorous mathematical 
well-posedness theoretical ground yet for these kind of boundary conditions, as 
opposed to the simpler case of maximal nonnegative boundary conditions. 
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Appendix A. Maximal nonnegative boundary conditions for 

SYMMETRIC HYPERBOLIC SYSTEMS 

Let n C M" be a bounded domain with smooth boundary and T > 0. We 
introduce the notations O = (0, T) x and F = (0, T) x dft. Consider the first order 
differential operator L := dt + J^^^i Ai{t, x)di + B{t,x), where Ai eLip(O), B e 
i°°(0), and {B + B*)/2~ J2i=i diAi £ L°°{0). We suppose that L is symmetric, 
that is, Ai = A* on O. Our interest is in solving the initial-boundary value problem 

Lm = /(t,.T) in O, u(0,-) = g in r2, u{t,x) e N{t,x) for {t,x) e [0,T] x dfl. 

(40) 

where N{t, x) is a Lipshitz continuous map fi-om [0, T] x dft to the subspaces of 
C". Set be the outer unit normal to T, and denote by An{t,x) the boundary 
matrix/operator An{t, x) := X]r=i We assume that F is characteristic 

of constant multiplicity in the sense that dim ker An is constant on each component 
of F. We next suppose that N is maximal nonnegative, that is, the following two 
conditions hold on F: 

(A„(t, x)v, v) > 0, V(t, x) e F, Vw e N{t, x) (41) 

and 

dim N{t,x) ~ nonnegative eigenvalues of An(t,x) counting multiplicity. (42) 

The maximality condition (|42p implies that the boundary subspace N cannot be 
enlarged while preserving (|4ip . 

Let Ho := {u G ^^(O) : Lu e L'^{Q)\. It is easy to prove that Ho is a Hilbert 
space with respect to the inner product (u,it)H,, :— {u,u)]j2((y^ + {Lu,Lu) jjnp,-^. 

Theorem 10. (L^ well-posedness. Theorem 9 in [51]) For any f G L'^{{0,T) : 
L^(f2)) and g G L^(r2) there is a unique u G Ho satisfying (j40p . /tt. addition, 
u£C{{G,T):L\n)), 

sup ||u(t)||L2(f2) < C'||/|jii((o,T):L2(n)) + ||5lU2(n), 

0<t<T 

and 

\\u{t2)\\mn) - \\u{ti)\\L2(n) < / \\f{cr)\\L^n) + C||7i(cr)||L2(n) da. 
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